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On the basis of ab-initio simulations, the value of strength of interatomic bonds in one-, two- 
and three-dimensional carbon crystals is obtained. It is shown that decreasing in dimensionality 
of crystal gives rise to nearly linear increase in strength of atomic bonds. It is ascertained that 
growth of strength of the crystal with a decrease in it dimensionality is due to both a reduction 
in coordination number of atom and increase in the angle between the directions of atomic bonds. 
Based on these data, it is substantiated that the one-dimensional crystals have maximum strength, 
and strength of carbyne is the absolute upper limit of strength of materials 
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I. BACKGROUND 

According to the existence paradigm of strength, there 
are two ways to increase strength, namely: (i) by increase 
in the density of defects in crystal lattice (increase in a 
lattice distortion), or, (ii) vice versa, creation of defect- 
free crystals. The first way is a typical approach to pro¬ 
duction ordinary “bulk” materials. High levels of the 
strength of nano-sized crystals are due, primarily, to the 
absence of defects — . However, recently it was shown 
that extremely high levels of the strength and modulus 
of elasticity can be reached by decrease in dimensional¬ 
ity of crystals i.e., by the transition from three- to two- 
and one-dimensional crystals — . Most clearly this ef¬ 
fect can be demonstrated for allotropic forms of carbon. 
Thus, the strength of defect-free 3-D placeCityCrystal 
(diamond) is 95 GPa — , of 2D-crystal (graphene) - 130 
GPa— and of ID-crystal (carbyne) - more than 270 GPa 
— . In this connection, the task of the work was to es¬ 
tablish the regularities of influence of the crystal dimen¬ 
sionality on it strength and to develop ideas about the 
physical nature of this effect. 


II. METHODS 

One-, two- and three-dimensional carbon crystals were 
used as an object of investigation. To determine the 
strength of a one-dimensional crystal (monatomic chains 
of carbon atoms), ab-initio calculations were employed. 
In this case, tension of the chain was simulated. Strength 
was determined for polyyne structure, which is energy- 
wise more favorable. Full energies of chains were calcu¬ 
lated by the pseudo-potential method (software package 
Quantum-ESPRESSO (QE)) — . This method was used 
for modeling the mechanical properties of the chains. 


Pseudo-potentials for carbon were generated under the 
scheme Vanderbilt Ultrasoft with the package Vanderbilt 
eode, version 7.3.4—. To determine the accuracy of the 
calculation of the total (full) energies of carbon chains, 
test calculations for infinite chains with the structure of 
polyynes were executed. Interatomic distances and total 
energies, consistent with the results of the work — were 
obtained. The value of the cutoff energy Eeut = 450 
eV. In the calculations, the exchange-correlation poten¬ 
tial PBE — was used. 

To determine the strength of atomic bonds, R 2 D 1 
two-dimensional (2D) crystal (graphene), the value of 
fracture stress of graphene — obtained by ab-initio sim¬ 
ulation of tension of graphene sheets in two directions 
“zigzag” and “armchair” were used. The magnitude of 
atomic bond strength Ee was determined as the value of 
critical force acting on the bond at the time of instability 
of the graphene sheet at it tension: 

R 2 D = <^dobo X (I - 5 ) (I) 


R^f^ = a'^'^'^dobo x (I — 5 ) x (I + cos (Oq)) x cos (I — 6>') 

( 2 ) 

where and are the critical stress of insta¬ 

bility of graphene sheet under tension in the directions 
’’zigzag” and ’’armchair”, respectively; bo = ^/3/2ao; 
oq = 0.1422 nm is the lattice parameter; do = 0334 nm 
is the effective thickness of graphene sheet; Oo = 120° 
is the angle between atomic bonds ij and jk in unde¬ 
formed state; 0' = 130° is the angle, at which instability 
in graphene occurs (was determined based on the geom¬ 
etry of the system at instability); 5 is the transverse de¬ 
formation of graphene sheet. 

The strength of atomic bonds in three-dimensional 
(3D) crystal (diamond) was determined on the basis of 
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ab-initio simulations of tension of diamond crystals in 
the direction < 111 > — . The strength value was cal¬ 
culated as the critical magnitude of the force acting on 
the interatomic bond at instability of the crystal under 
tension: 


Rsd = ^^CTcalcos"^ (7r/2 - 9) (3) 

where dc is the critical stress of instability , oq is the 
lattice parameter, 0 is the angle between atomic bonds 
ij and jk in deformed state. 

Theoretical analysis of the orientation dependence of 
strength and its value dependence on coordination num¬ 
ber of atom was executed within the formalism of Bren¬ 
ner potential—, which is based on Abel pseudo-potential 
theory — that describes well the carbon compounds. Pa¬ 
rameterization of the potential was carried out by the set 
II of Brenner parameters—: R = 0.139 nm, D{e) = 6.0, 
5 = 1.22, 13 = 0.21 (5 = 0.5, uq = 0.00020813, 

Cq = 330^, dg = 3.5^. Since the Brenner potential in ini¬ 
tial state gives a significant overestimation of strength of 
atomic bonds, which reaches values of 30nN (Fig. [TJi) , 
its modification was carried out for the crystal that con¬ 
sisted in correction of cutoff function fcuti^ij) = 1. It 
enabled to get the values of strength, which agree well 
with the results of ab-initio simulations. 


III. RESULTS AND DISCUSSION 

The strength of atomic bonds in one-dimensional crys¬ 
tal was determined by the critical stress of instability of 
monatomic chain under the conditions of uniaxial ten¬ 
sion (Fig. [2]). In accordance with the data obtained. 
Rid = 11-3 nN. This value is in good agreement with 
the evidence of 12.2 nN — , and with the value 12.3 nN, 
obtained for chains containing more than ten atoms — . 

To determine the strength of atomic bonds in 2D- 
crystals, data of ab-initio simulation of tension and frac¬ 
ture of graphene sheets in ’’zigzag” and ’’armchair” di¬ 
rection — were used. The critical value of the force at 
the moment of instability of the graphene sheet was de¬ 
termined by equations (HD and ©• According to the re¬ 
sults obtained, the strength of atomic bonds in graphene 
for ’’armchair” direction was 8.9 nN and for ’’zigzag” di¬ 
rection - 8.3 nN. It should be noted that the difference 
in angles between atomic bonds in the stretched sheet 
of graphene is the cause for different values of atomic 
bonds strength in graphene at transition from ” zigzag” to 
’’armchair” orientation. Thus, at the moment of instabil¬ 
ity under tension in ” armchair” direction the orientation 
of atomic bonds is characterized by angles of 130° and 
130°, while for tension in ’’zigzag” direction, they are 
114°, 132°, respectively (Fig. [3)). Really, calculations 
using Brenner potential, in the first case give the value 
of strength of atomic bonds equal to 8.49 nN, and in the 
second case this value is 7.5 nN (Fig.[3l). Accordingly, the 



FIG. 1. (Color online) The original and modified Brenner 
potential. Force vs. strain (). Energy vs. strain (b). 



FIG. 2. (Color online) Dependence ’’force vs. strain” for 
infinite monatomic chain. Result of ab-initio calculations. 
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FIG. 3. (Color online) Dependence of the strength of inter¬ 
atomic bonds of carbon crystals on the its dimensionality: ■ 
are the strength values obtained by the results of ah-initio 
calculations; and • are the results of calculation by the mod¬ 
ified Brenner potential for angles between interatomic bonds 
in the unstrained state and at the critical strain of crystal 
instability, respectively. 


ratio of the strengths of atomic bonds obtained accord¬ 
ing to ah-initio calculations is 1.07, and of those obtained 
from the Brenner potential is 1.14. The difference does 
not exceed 9%. The absolute values of these strength 
’’anisotropy” are not much great, but they illustrate one 
of the specific features of interatomic interaction of car¬ 
bon atoms, namely, its orientation dependence. 

The value of bond strength in diamond was obtained 
based on the results of ab-initio simulations. The criti¬ 
cal value of the force acting in the bonding at the time 
of instability of the diamond lattice was determined by 
the formula ([3j). In compliance with the results obtained 
bond strength in diamond is Rsd = 5.2 nN. 

According to the results, change in the dimensional¬ 
ity of defect-free crystals essentially effect on their bond 
strength. This effect manifests itself in a significant 
(more than 2 times) increase in the bond strength of 
the transition from three- to two- and one-dimensional 
crystals (Fig. [3]). Strength values of bonds calculated 
using the modified Brenner potential are plotted on the 
same graph. Calculations were executed for both the 
initial crystal and taking into account changes of angles 
0 between the bonds when reaching a critical strain of 
crystal instability. The results of calculations, agree well 
with the ah-inito results. As it is indicated in Fig. [H 
the transition from three-dimensional to two- and one¬ 
dimensional crystals, change in coordination numbers is 
accompanied by a change in the angles 0 between the di¬ 
rections of atomic bonds. To separate the contributions 
of the coordination number Z and the angle 0 to this 
change in the strength, i?, Brenner formalism was used. 
In explicit form the dependence of strength R on the co¬ 



FIG. 4. (Golor online) Influence of bonds orientation on the 
strength of carbon crystals at a fixed value of coordination 
number Z: • angles between the bonds at the moment of 
instability of carbyne (180°), graphene (130°) and diamond 
(113°). 


ordination number 2 : may be derived only for the case of 
equal angles Oijk with all neighboring atoms k: 

R{Z, 9) = - sT^) (1 + (^ - 1)G (Oijk)) 

( 4 ) 

where Gi{0ijk) is angle function—. 

Using the relations suggested in — and the magnitudes 
of Brenner parameters (set II from—), orientation depen¬ 
dencies of strength at fixed values of coordination number 
Z = 2, 3,4 were built (Fig. |4]). 

Formally, Brenner dependence describes the change in 
strength with the angle within the range from 0° to 180°. 
But in the real objects (diamond, graphene, carbyne) 
variation of the angle does not exceed the range 90° —180° 
deg, so the analysis is confined to this range of values. 

According to the data shown in Fig. 01 the sensitiv¬ 
ity of strength to change in the angle increases with the 
growth of coordination number Z. On the other hand, the 
sensitivity of strength to changes in coordination number 
increases with decreasing in angle. As it is indicated in 
Fig. [5l at transition from diamond to carbyne, bond 
strength increases 2,3 times, while 56% of this increase 
in the strength is due to a decrease of the coordination 
number and 44% is related to the growth of angle to the 
maximum possible one, which is 180°. 

Regularities of dependence of atomic bond strength 
on coordination number is demonstrated by the carbon, 
because it is most convenient to analyze this problem. 
However, these regularities are inherent also to the other 
materials. Work — gives evaluation of the strength of 
atomic bonds in ID and 3D crystals for metals such as 
Cu, Ag, An, Pd and Pt. The results of calculations by 
different methods shown that at transition from ID- to 
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FIG. 5. (Color online) Dependence of the atomic bond 
strength on the coordination number at the fixed values of 
angles. A are the critical value of angles. 


3D-crystals, the strength of atomic bonds should increase 
by 2 — 3.3 times. Thus, the material in the form of a one¬ 
dimensional crystal (monatomic chain) should have the 
highest level of strength. In this regard, strength of car- 
byne claims to be the absolute upper limit of strength 
of materials. According to the results of ab-initio cal¬ 
culations, carbyne containing 5 atoms has a maximum 
strength. This strength is 13.1 nN — . When the effec¬ 
tive diameter of the chain equals 0.200 nm—, its value is 
417 GPa. In — for an infinite chain the strength values 
9.3 — 11.7 nN were obtained, and considerably smaller 
value of effective diameter d = 0.077 nm was utilized. 
This gives the strength value for infinite chain equal to 
1988 — 2501 GPa. In — the experimental measurement of 
carbyne strength was executed. It was ascertained that 
its strength must exceed 270 GPa. This lower bound 
for strength of carbyne is more than 2 times greater 
than the experimental value of strength of graphene (two- 
dimensional crystal of carbon). 


IV. CONCLUSIONS: 

Increase in strength of the interatomic bonds while re¬ 
ducing in dimension of the crystal is a specific feature of 
nano-sized defect-free crystals. At transition from three 
to two and one-dimensional crystal of carbon, interatomic 
bond strength increases almost linearly. The bonding 
strength in carbon monatomic chain is 2.2 times higher 


than that in diamond. This value is of the same order as 
the increase in strength of metallic crystals, which may 
be 2-3 times. 

Growth of crystal strength with decreasing in its di¬ 
mensionality is due to two factors, namely, a decrease in 
coordination number (the number of atomic bonds) and 
increase in the angle between the bonds. At transition 
from diamond to carbyne, 56% of increase in strength is 
due to a decrease in coordination number and 44% - due 
to increasing in angle between the bonds. 
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